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I attempt to quantify how far from maximal one should expect the atmospheric mixing angle to be 
given a neutrino mass-matrix that leads, at zeroth order, to a vz mass-eigenstate that is 0% v e , 50% 
Vp, and 50% v r . This is done by assuming that the solar mass-squared difference is induced by an 
"anarchical" first order perturbation, an approach than can naturally lead to experimentally allowed 
values for all oscillation parameters. In particular, both | cos 2# a tm | (the measure for the deviation of 
atmospheric mixing from maximal) and \Uez\ are of order -j Am^/Am^ in the case of a normal 
neutrino mass-hierarchy or of order Am s 2 ol /Am^ tm in the case of an inverted one. Hence, if any of 
the textures analyzed here has anything to do with reality, next-generation neutrino experiments 
can see a nonzero cos 2# atm in the case of a normal mass-hierarchy, while in the case of an inverted 
mass- hierarchy only neutrino factories should be able to see a deviation of sin 2 26 atm from 1. 

(N ; 

X> 1 I. INTRODUCTION 

: 

One of the most exciting developments in particle physics of the past six years has been the confirmation, by several 
different experiments 0, 0i 13 ■ that neutrinos have mass and that, similar to quarks, neutrino mass-eigenstates and 
weak-eigenstates mix. The existence of neutrino masses is, presently, the only solid, unambiguous evidence of physics 
1 that cannot be accommodated by the standard model of electroweak and strong interactions. However, while there 
^ ' is evidence for physics beyond the standard model, the amount of information available is very scarce. The two most 
significant facts we did learn (and would like to understand better) are: i-neutrino masses are significantly smaller 
than all other fermion masses (m u /m e < 1CP 6 ), ii-the leptonic mixing matrix is qualitatively different from the quark 
mixing matrix. In this paper, I'll have nothing to add regarding (i), but would like to address (ii). 

The quark mixing matrix is well approximated by the unit matrix (see, for example, Q), i.e., the mixing angles of 
the CKM matrix are either small (sin 2 9c — 0.05) or very small (| V^l 2 ~ 10 -5 ). Combined with the fact that quark 
. masses are hierarchical (nit 3> m c 3> m u , etc), the peculiar form of the CKM matrix has led to the belief that there 
is some underlying symmetry or dynamical mechanism that fundamentally distinguishes the different quark families, 
and naturally explains the small mixing angles and the disparate mass eigenvalues. Indeed, it seems quite implausible 
that the CKM matrix is representative of a "random" unitary mixing matrix.* 
Oh, In the leptonic sector, the situation is quite distinct. The leptonic mixing matrix is "far" from a unit matrix, and 
all its entries are large, with the exception of the U e 3 element,^ which is somewhat small: |J7 e 3| 2 < 0.07 at three sigma 
confidence level 0,0. It is not clear that there is, currently, strong evidence to believe that a new underlying symmetry 
or dynamical mechanism that distinguishes the different leptonic families is required to explain the leptonic mixing 
angles. 1 '" Indeed, it has been proposed that perhaps the neutrino-mass sector of the standard model is structureless, 
and that the leptonic mixing matrix is a reasonable representative of a random mixing matrix |10| . 

This statement has recently been verified in a more quantitative way |5j. It was pointed out that the "anarchical 
hypothesis" for the leptonic mixing matrix is consistent as long as \U e s\ does not differ significantly from its current 
upper bound. Furthermore, if correct, the anarchical hypothesis provides a lower bound for |J7 e 3|> namely |t/ e 3| 2 > 
0.011 at the two sigma confidence level This means that next generation searches for |t/ e 3| 2 using, say, neutrino 
"super-beams" ^l) or reactor antineutrinos ^1 can definitely test the hypothesis that the leptonic mixing matrix is 
consistent with a random unitary matrix. 

In spite of the claims above, a significant amount of research effort has gone into determining what sort of underlying 
symmetries or mechanisms can explain the leptonic mixing angles and the neutrino masses jl-ill . There are several 
reasons for this, and I will spell out two of them. One is the belief that, at some very high energy scale, the distinction 
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* An attempt to quantify this statement can be found in 

t I define the mass eigenstates such that m| > m 2 , and Am 2 3 23 1 > Am 2 2 , where Am 2 . = m 2 — m? . Weak-eigenstates u a are related 
to the mass-eigenstates Vi via v a = U a iVi, where a = e, fi,T , i = 1,2,3, and U a i are the elements of the leptonic mixing matrix, also 
referred to as the MNS matrix, Umns- Mixing angles 9ij are defined as per the PDG prescription Q, 

t One potential criticism is that the charged lepton masses are hierarchical. While this issue will not be addressed here at all, it is easy 
to imagine scenarios where the charged lepton masses are hierarchical, while the leptonic mixing matrix is structureless Indeed, in 
some scenarios it is possible to obtain a structureless mixing matrix while charged lepton and neutrino masses are hierarchical • 
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between quarks and leptons is erased. This happens, for example, in most grand unified models. Hence, because 
the quark sector seems to be structured, it is likely that the same is true of the leptonic sector [lj] (for a counter- 
argument, see @)- Another reason for believing in structure in the leptonic section is the curious fact that, currently, 
the atmospheric mixing angle is consistent with maximal (823 — 7r/4). While maximal atmospheric mixing is in 
agreement with the anarchical hypothesis |5j, it can also be interpreted as evidence for specific "textures" in the 
neutrino mixing matrix, and most of these textures can be related to underlying flavor symmetries. 

In this paper, I will try to address how close to maximal the atmospheric mixing angle must be if its "maximality" 
is indeed related to some fundamental principle. In the next section I discuss what assumptions will be made, and 
define how this issue will be addressed. After that, I'll analyze different mass-matrix textures that lead to maximal 
atmospheric mixing at leading order, and estimate the typical deviation of the atmospheric angle from maximal, 
among other observables. Discussions, a summary, and conclusions follow. 



II. NEUTRINO MASS MATRICES, MAXIMAL MIXING, AND RANDOM PERTURBATIONS 



Instead of concentrating on several distinct flavor models and dynamical mechanisms for explaining the pattern 
of neutrino masses and the values of the mixing angles, I'll concentrate on a bottom-up approach which imposes, 
at the leading order, maximal mixing in the atmospheric sector. For this purpose, I select zeroth order neutrino 
mass-matrices Mo that, in the basis where the charged current and the charged lepton mass matrix are diagonal, have 
one mass-eigenstate which is 0% v e , 50% and 50% v T . I will only consider Majorana mass matrices and do not 
particularly care about the "fundamental origin" of Mq. The next step is to perturb Mq by adding a generic "first 
order perturbation" mass matrix Mi. See the Appendix for a definition of "generic." As in the case of Mo, the origin 
of Mi is not relevant to the considerations made here. It can be associated to the renormalization group running of 
the several couplings in specific models and/or heavier see-saw right-handed neutrinos and/or a different, sub- leading, 
source of lepton-number violation, etc. 

In order to make more quantitative assertions, it is important to define the size of the perturbation. This will 
be done by requiring that, after the perturbation is added, other features of neutrino mixing (especially the solar 
mass-squared difference) are properly explained. 

It proves instructive to start with two toy 2x2 scenarios (for concreteness, the x v T "subspace"). The first one 

is 
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where mo, mi (too S> mi) have dimensions of mass and a,b,c are O(l) real numbers. § It is easy to see that the 
eigenvalues of Eq. (|II.1(I are mo,— mo, and its eigenvectors are (l/v^, l/^/2) T and (l/v^, — ^/V%) T , respectively. 
While it does not concern the discussions here, it is worthwhile to mention that Eq. (|II.1(I can be interpreted as due 
to an — L T global symmetry. 

It is trivial to diagonalize Eq. (|IL1J| +Ea. I|II.2J) exactly and obtain eigenvalues and eigenvectors. It proves useful, 
however, to obtain approximate results via perturbation theory. Explicitly, 
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(II.4) 
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Here, a, (5 = ^(a + b T 2c) and e = ^.(o - b). From Eq. 



1/V2 1/V2 J 
9 it is easy to read off the mass-squared splitting 



§ It suffices to consider only real mass matrices, which is assumed henceforth. One can always deal with complex mass-matrices by 
diagonalizing MM* = C/diagonal[|mi | 2 , \rr12 1 2 , • • ■]U* instead of M = C/diagonal[mi, mi ■ ■ -]U T , and obtain, qualitatively, the same 
results. 
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induced by the perturbation Mi, 

Am 2 



(-1 + a) 2 - (1 + (3) 2 \ ~ 2 |a + /3| = 2^ |a + b| , (II.6) 
1 m 



"'0 

(higher order terms have been dropped) and the leptonic mixing matrix, 

"-'-^(-r-'A-i)- < iL7 » 

One way to quantify how far from maximal is the mixing angle is to compute | cos 20 \ = | cos 2 — sin 2 0\. From the 
expression above, ignoring higher order terms, 

I «*l - 5 Id - «>' - & + >?\ = * 14 - ^ l« " 61 = (IL8 » 

What do these results mean? If these are compared to the experimental data, it is natural to choose m 2 , to be the 
atmospheric mass-s quared, Am 2 tm ~ 2 x 10~ 3 eV 2 (see \n£ for a recent three-flavor analysis of the combined solar, 
atmospheric, and KamLAND data) and to associate the induced Am 2 , Eq. (|II.6|1 to Am 2 ol ~ 7 x 10 -5 eV 2 , and to 
the atmospheric angle sin 2 (9 2 3 ~ 0.5. Hence, assuming \a + b\, \a — b\ both to be order one, 

m ~ 4 x 1CT 2 , (II.9) 
mi ~ 9 x 10~ 4 , (11.10) 
| COS 2023 1 <0.01. (11.11) 

Next, consider the so-called "democratic" texture, namely 

mo = ^ ( ; ; ) , (H.12) 

with eigenvalues m ,0 and eigenvectors (l/v 7 ^, 1/V%) T and (1/V%, — 1/V%) T , respectively. In this case, Eq. I)II.12|1 + 
Eq. 1II.2|) is, ignoring higher order terms (this is very similar to the previous case, Eq. Jn3Jl), 

where a, (3 = ^^-(a + & T 2c), e = (a — b). The interpretation of the different terms in this case is somewhat 
different. As before, one can naturally associate m 2 , to Am 2 tm while it is natural to associate the shift of the zero 



eigenvalue of Eq. (|II.12(I to a/A 



sol * 



Am 2 tm 4m 2 



1 I (a + 6-2c) 2 . (11.14) 



The modified atmospheric angle is, hence, 



mi , , / Am 2 , 2 la — 61 , , 

cos26» 23 = —\a-b\~ \ - 9 SQl -. — ^ i— <0.2, 11.15 

' ' rn l 1 V Am 2 tm |a + 6-2c| ~ V j 

where, in the last step, |a — b\,(a + b — 2c) 2 where replaced by 1 and 37r/2, respectively, their average values.^ 

Eqs. Ijll.lllll,15|l will be considered as the upper bounds on | cos 2023 1 for the two different "models." Note that 
while the two Mq yielded the same zeroth order physics, Eq. (|II.11|) and Eq. (|II.15(1 differ by an order of magnitude! 
The reason for this is simple: in the case of Eq. I|II.15|) . the perturbation parameter mo/mi was chosen much larger 
than in the case of Eq. (|II.11(I in order to fulfill the same intuitive agreement with the "sub-leading" oscillation 
parameters. It will become clear in the next section that this discussion of 2 x 2 toy models does indeed capture the 
relevant features of more realistic 3x3 mass matrices. 



■ See the Appendix for the prescription for combining order one "random" factors. 
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III. CASE STUDIES 

In this section, several examples of 3 x 3 mass-matrices Mq + Mx that yield, at zeroth order, maximal atmospheric 
mixing are inspected. Some of these have already been identified in previous studies, and have been used to determine, 
for example, expected values of \U e3 \ for different neutrino mass textures 0|- Different M differ by the values of 
other parameters, both known and unknown, such as the solar mixing angle and the hierarchy of the neutrino masses. 
The approach used in the previous section will be employed in order to estimate the maximal | cos 2# 2 3 1 that each M Q 
allows given a random perturbation Mi. The philosophy of studying the effect of random, sub- leading perturbations 
on flavor structures was already followed, for example, in [lfj. 

Before proceeding, it is worthwhile to comment on | cos 2023 1 &s the relevant measure of deviation- from-maximal. 
Maximal atmospheric mixing means |t^3| = 1/2, such that deviation from maximal is proportional to | 1 2 — 1/2 = 
— 1/2 cos 2#23 — sin 2 #13 sin 2 #23 ■ Note that the second term in the previous expression is proportional to |J7 e 3| 2 , much 
smaller than the first term (this can be easily verified a posteriori for all mass-matrices discussed here), such that 
I 1 2 — 1/2 ~ —1/2 cos 2^23- Hence, cos 26*23 is a bona fide measure of deviation- from-maximal-atmospheric-mixing. 

Analogous to the previous section (Eq. i|II.2|n . 




M x = m x \ d b f , (IIL1) 
\e I cj 

will be assumed for all the discussions that follow, a, b, c, d, e, / are random order one real numbers. 

A. Normal Hierarchy 

A "normal" mass- hierarchy is characterized by Amf 3 > 0, where to 2 , 3> to 2 , to 2 .. There are several mass-matrix 
textures that lead to two massless and one massive eigenstate, the latter composed of a 50-50 — v r mixture. 
One possibility is to introduce one right-handed neutrino with a very heavy Majorana mass, coupled with identical 
strength to L M and L T * plus the Higgs doublet, in the basis where the charged lepton Yukawa couplings and the 
charged weak gauge couplings are diagonal |9l Il7|. Upon integrating out the right-handed fermion, and below the 
electroweak symmetry breaking scale, 
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l/y/2 -l/y/2 I (III.2) 



is generated. Defining 



-1/V2 1/V2 j \ 1 J \ 1/V2 1/V2 



1 

U = \ l/y/2 l/y/2 I , (III.3) 



-l/y/2 l/y/2 _ 



Mq + Mx can be written as 



a e 
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M + Mx = m Uo \ e r? | Uj , (HI.4) 
S 77 1 + 7 



where 



mo ' " mo V 2 J J * 1 mo V 2 J / 



e= rrn±_e g = mi d±e »H k_£ ( m 5) 

mo V2 ' mo V2 ' ' mo 2 V ' 



hi are the left-handed lepton weak doublets, (i^,^) 1 
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Therefore, at leading order in mi/mo, 



cos a; sinw S \ / a' 
Mo + Mi = m £/o ( — sin a; cosw 77 J I /3' ] ( sinw cosw —7/ ] C/q 1 ^. (III. 6) 

-5' -7]' 1/ \ 1 + 7 
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a' and /?' are the eigenvalues (ordered in ascending order of absolute value) and u> the mixing angle of the 2x2 sub- 
a e 
e 



matrix ( a 1 I , while 5' = S cos w — 77 sin w and rf = 5 sin lu + rj cos w. The two independent mass-squared differences 



arc: 
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while the mixing matrix is 



cos a; sin a; <5 

tfatfs = I * * 1/V2(1 + »»)), (nii 



1/V2(1 




where only the relevant elements have been spelled out (others have been mercifully represented by *). One can then 
"read off" 

012 =w, \u e3 \ = \s\ = ^^, (IH.9) 

I cos20 23 | = §| (1 + v f - (1 - - 2 M = ^ l b - c l ■ ( IIL1 °) 

As done in the previous section, estimates for the sizes of various parameters are obtained after mi /mo is extracted 
from the fact that the ratio of mass-squared differences is known: Am 2 2 /Am 2 3 ~ 0.034. More quantitatively, it 
has been estimated that, at the three sigma confidence level, 0.018 < Amj 2 /Amf 3 < 0.054 [l5|. A revised analyses 
of the Super-Kamiokande at mosp heric data, however, seems to point to a slightly smaller value for the atmospheric 
mass-squared difference (see |l8j . There is, currently, no publication of this reanalysis). In order to take this into 
account, I simply rescale the results of by the ratio of the "old" and "new" central values (2.6 x 10 _3 /2.0 x 10~ 3 ): 



0.023 < r- m * 2 , < 0.070, (III.ll) 

A™ 2 1 
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which should be very close to the real three sigma confidence level range. For all estimates bellow, in order to be 
conservative, I'll use the upper bound of Eq. (|lll.ll[l . which yields the loosest constraints on cos2# 23 . In the case at 
hand 



'12 



_ \j3 + a \ V(/3 - a) 2 + 4e 2 = -\ 
Ato 2 3 1 1 v y ' m\ 




(III. 12) 



Replacing the quantity inside the absolute value by its average value (see Appendix) , 

to? 0.070 



1.5 



(0.2y, (111.13) 



such that 



\U e3 \ < -0-09, |cos20 23 |<O.2W|~O.18. (111.14) 



t The PDG definition of the leptonic mixing matrix is u a = U a pVi for a = e,[i, r and i = 1,2,3, and the mass matrix in the flavor 
basis is given by U* diagonal [mi, W2 , ms]U'. Hence, Eq. 1111.81 is equal to U*. However, since I am restricted to real mass matrices, 
U = U* and this distinction can be safely ignored. 
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In summary, it is expected that the solar angle is "anarchical," i.e., dictated by the diagonalization of a 2 x 2 sub- 
matrix with random order one entries. This hypothesis fits the data very well 5j. On the other hand, both \U e s\ and 
the deviation of 6*23 from maximal are of order a/Ato^/Ato 2 ^, O(10 _1 ). 

A couple of important comments are warranted. The average upper bound on | t/ e 3 1 2 < 0.008 is smaller than the 
estimated lower bound on |t/e 3 | 2 obtained from the anarchical hypothesis |5J (albeit only slightly), meaning that 
next-generation experiments sensitive to \U e 3\ 2 ^ 10~ 2 may be able to start distinguishing the two paradigms. The 
average deviation of 823 from maximal is rather large, but well within current bounds |15| : 

I cos26> 23 | < 0.34 (0.44), at 2a (3a). (111.15) 

Finally, it should be noted that given the approach adopted here, nothing can be said about the sign of cos 26*23, 
i.e., whether the 1/3 state is "predominantly (cos 26*23 < 0) or "predominantly v T " (00326*23 > 0). Note that we 
already know, experimentally, that V3 is not "predominantly v e " (\U e 3\ 2 < 0.5) and that, thanks to matter effects in 
the Sun's core, V\ is "predominantly v e " (cos 26^2 > 0, the "light-side" of the solar neutrino parameter space ^{|) at 
the five sigma confidence level |l5j . 



B. Inverted Hierarchy 

An "inverted" mass-hierarchy is characterized by Ato 2 3 < 0, and two eigenstates that are almost degenerate in 
mass-squared. As in the normal hierarchy case, there are several textures that yield two massive mass-eigenstates, 
degenerate in mass-squared, and a massless one, composed of a 50-50 v„ — v T mixture. I'll look at three different 
cases, characterized by whether the massive mass-eigenstates are degenerate in mass (as opposed to mass-squared) or 
by how the massive mass-eigenstates "share" the v e content. 

The case which is "most similar" to the one in the previous subsection is well represented by 
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-1/2 1/2 

Note that the heavy mass-eigenstates are degenerate in mass. Mq + Mi can be written as 



UJ . (111.16) 




M + Mi = m U e 1 + j3 r, t/ ' , (111.17) 



where a, (3, 7, S, e, S, and n are as in Eq. I|III.5(I . 

As in the previous case, degenerate perturbation theory yields 

(cosw sinw — S \ / 1 + a' 0\ / cosw — sinw 5' \ 
-sinw cosw 1 + /?' ] sinw cosw rf Uj , (III. 18) 

5' i 1 J \ 7 J \ -5 -T) 1 J 

where the primed quantities and lo are defined as in the previous subsection. 

The extraction of oscillation parameters proceeds as before. Indeed, the MNS matrix is identical to Eq. (jIII.8|l . up 
to 5, r] — * —<5, —rj, while the mass-squared differences are 

Am 2 

Am 2 3 = -m 2 , 1^ = 2 (/?'-«')• (111.19) 

Note that Amfg also receives a higher order (O(momi)) correction, which can be safely neglected. Furthermore, As 
before, we extract mi /too by requiring 



^ = 2^(P~af+4e2 = 2^J(a- b -±^ + f) + 2(d - e)» * 0.07. ,TTT.20) 

ATO^g " TOQ V \ & 



where I conservatively use the upper bound of Eq. (fill - 1 1 1) . One obtains, on average, mi/mo — 0.02, which implies 

02 [2 [2 
l^eal <— Yg - 0.008, I 00820231 < 0.02^- -0.016. (111.21) 

Hence, in the case of an inverted mass-hierarchy with two mass-degenerate (at leading order) eigenstates, one expects 
a large, "random" solar angle, and very small values for |£7 e3 | and cos 26*23, both of order Ato 2 2 /|Ato 2 3 | (O(10 -2 )), 
well within current experimental bounds. 
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C. Inverted Hierarchy, Traceless 

There is the possibility that the two massive zeroth-order mass-eigenstates have opposite CP-parities (i.e., mo,i = 
-771(^2 )• One example is 

/ 1 \ / 1 \ 

T 











j = m U ( 


-1 












M = m -1/2 1/2 =m U 0-10 U '. (111.22) 



Similar to the previous two cases, Mq + Mi can be written as 




..!/„ + V/, = m.J-n I 6 -1 I )' // 1 r (l - (111.2:!) 



where a, /3, 7, (5, e, S, and 77 are as in Eq. (|III.5|) . 

In this case, the diagonalization procedure is somewhat different, as there are no degenerate eigenvalues (it is, 
indeed, simpler). To lowest order in mi/mo, 

1 -e/2 -S\ f 1 + a 
M +M 1 = m U I e/2 1 77 -1 + /3 | | -e/2 1 -77 | {/„' . (111.24) 

S -77 1 J \ 7 

The two independent mass-squared differences are, to lowest order: 

Aml 3 = -ml i^ = 2|/3 + a| = ^|2a+& + c-2/|, (111.25) 
|Am 13 | mo 




while the mixing matrix is 



Hence 



-e/2 -5 

Umns = ( * * 1/V2(1 +»?)). (HI.26) 

* 1/72(1-7?) 



! sin ^l = I = ^t-J d - e l> l^sl = 1*1 = (HI.27) 
|cos20 23 | = i|(l + 7 7 ) 2 -(l-7? 2 )| =2|r ? | = ^|&-c|. (111.28) 



Following the by-now-familiar procedure, one estimates 

. , Am?, V2 , Am?, 1 , , Am?, 1 

l a '°""l S iA^ivf I^ISjAsfiVS' i^^iSjAsfivi- <'"- 29 > 

Note that, here, not only are \U e3 \ and cos20 2 3 smaller than Am 2 2 /|Am 2 3 |, but so is the solar angle (indeed, it is, 
on average, smaller than #13!). Given the current lower bound on the solar mixing angle, sin 2 # 12 > 0.22 at the three 
sigma confidence level MTone can only conclude that such a scenario is experimentally ruled out, as the estimate 

above is sin 2 6 12 < 0.0005. 

One should not, however, conclude that the mass-matrix texture Eq. 1111.22(1 is currently ruled out, but that if it is 
to work, it has to be "perturbed" by a rather "structured" matrix in order to guarantee that the solar angle is very 
strongly modified from its leading order value, while, say, |J7 e 3| is not. This could be accomplished, for example, if d 
is close to (— e) and much larger than a, b, c, f } 



* This should be done with care. Mi stops being a perturbation as d approaches rao/rai, indicating i-that the approximate results 
obtained here need not be reliable, ii-that the zeroth-order mass-matrix texture has to be revised in or der to somehow include the 
information regarding the large-but-not-maximal solar angle. Such examples can be found, for example in and references therein. 
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D. Inverted Hierarchy, Traceless, Bi-Maximal Mixing 



An interesting mass-matrix texture is motivated by an L e — — L T discrete symmetry. It is 

/0 1 l\ /1/V2 l/y/2 \ (I 0\ /1/V2 1/2 1/2 

M =-p 100 =mo 1/2 -1/2 -l/y/2 1 1/^2 -1/2 -1/2 
V 2 \ 1 / \ 1/2 -1/2 l/y/2 ) V / V -l/y/2 l/y/2 



It proves economical to define 

/ l/y/2 l/y/2 
= 1/2 -1/2 -1/V2 ) . (HI.31) 
V 1/2 -1/2 l/y/2 

Note that not only do the massive mass-eigenstates have opposite CP-parity, but they also "share" the v e content 
equally {\U%$\ 2 = \U^ 2 \ 2 = 1/2). Eq. (|III.31I) is often referred to as a bi-maximal mixing matrix [2l| (sin 2 6 23 = 
sin 2 612 — 1/2, sin 2 8i 3 = 0), and has received a significant amount of attention in the literature. 

/ 1+a e S\ 

M + Mi = m C/ frm e - 1 + rj £/ bmT , (111.32) 
\ 5 1] 7 / 

where 

a -^(2 a + 6 + C + 2v / 2(rf + e ) + 2/), /3=^-(2a + & + C -2 v / 2(d + e ) + 2/), 7 = ^ (*±s _ /) f 

e =i^(2«- & -c-2/), 6=^(-V2(b-c)-2d+2e), ry = ^ (V2(6 - c) - 2d + 2e) . (111.33) 

Identical to the previous subsection, to lowest order in mi/mo, 

l/o h.l/i /// C/ fcm I e/2 1 77 I I L+ S n ) I </2 1 •// )?/,"" . (TTU4) 




and the two independent mass-squared differences arc, at leading order: 



Am 2 3 = -m 2 , =2|/3 + a| = ^-|2o + 6 + c + 2/|. (111.35) 



The leptonic mixing matrix is 



l/V2(l + e/2) 1/^(1 -e/2) l/y/i^-S) 
Umns=\ * * 1/V2(l + 1/V2fa + *)) I , (IH.36) 

* * 1/V2(1 - l/yffln + S)) 



such that 



|cos20 ia |=e=-^|2o-&-c-2/|, | cos20 23 | = V2|* + »,| = ^l| e -d|, \U e3 \ = 2* = £*-|6_ c |. (111.37) 
4ra mo V2 2m 

Solving for mi/mo in terms of the ratio of the mass-squared differences and replacing the order one expressions by 
their average values 

Similar to the previous subsection, this scenario runs afoul of the experimental data, because it predicts an incorrect 
value for the solar angle. It requires, on average, that solar mixing is too close to maximal: cos20 12 < 0.04, while 
current data indicate 0.22 < cos2# 12 < 0.56 at the three sigma level As in the case of the previous subsection, 
this result does not imply that the texture Eq. (|III.30(1 is ruled out, but simply that the "perturbation" has to be such 
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that the solar angle is significantly modified from its zeroth order value, while other parameters are less "perturbed." 
This could be accomplished, for example, if a is of order (— /), and much larger than b,c,d,e (but see footnote in 
subsection HUH ). 

There are other zeroth order bi-maximal textures that do not suffer from a too-large 9i2- The mass-matrix 

/ 1 \ / 1 \ 

M =m 1/2 1/2 =m C/ 6m 10 £7 6mT . (111.39) 
y 1/2 1/2 J \ J 

has two degenerate eigenvalues, indicating that the zeroth order solar angle is very sensitive to a generic perturbation, 
no matter how weak (this was observed, for example, in subsection HUB ) . In this case, one expects that both \U e 3\ 
and | cos2#23| are small, of order Am 2 2 /|ATO 2 3 |, while the solar angle is "anarchical." 

Finally, I'll conclude with a few words on textures that yield three mass-squared-degenerate eigenstates, some of 
which can be found in |l3j . In this case, it turns out that it is possible to obtain, for example, an anarchical solar 
angle while maintaining a generically small cos 26* 23 and \U e s\. The biggest challenge, one that cannot be overcome by 
a generic perturbation of the type Eq. Ijlll.ll) . is generating simultaneously the two distinct mass-squared differences. 
This is very simple to understand. The three eigenvalues which are degenerate at zeroth order will be modified to 
mi = mo(±l + (i), where Q are order mj/mo. Hence, Am 2 2 ~ 2(£ 2 — Ci)> Am 2 3 — 2(C3 — CO; which leads to 
Atijj 2 ~ Ato 2 3 , in conflict with the experimental data. Note that in the case of degenerate zeroth order eigenvalues it 
is quite unlikely that both mass-squared differences are induced by a single random perturbation - significantly less 
likely than in the original anarchy proposal |l0j . In an "anarchical model" that claims to explain both the mixing 
matrix and the mass-eigenvalues as typical fluctuations of a random mass-matrix (see the requirement on the 
eigenvalues of the mixing matrix is that they are "separated" by a factor y/\Amf 3 \/ Amf 2 ~ 5. In zeroth-order- 
degenerate models, the requirement imposed on the eigenvalues of the random perturbation-matrix is that their ratio 
is of order |Amf 3 |/ Am 2 2 ~ 30. 

As before, degenerate mass-matrix textures are certainly not ruled out, but they require a somewhat "structured 
perturbation" (see 22] for an example) in order to properly accommodate the neutrino data. In some cases, the 
lcptonic mixing matrix can still turn out to be completely anarchical. 



IV. DISCUSSION, SUMMARY, AND CONCLUSIONS 



The fact that the atmospheric mixing angle is currently consistent with maximal is often interpreted as evidence 
that there is a fundamental symmetry or dynamical reason behind the observed values of lepton masses and mixing. 
Indeed, several flavor models have been explored in the past that yield, at leading order, maximal atmospheric mixing 
and zero U e s. Sub-leading effects, however, will almost always lead to departure from these leading order predictions. 
These include spontaneous or explicit symmetry breaking effects and renormalization group running, and are often 
responsible for explaining the "solar" part of the neutrino oscillation parameters. 

It is important to understand how close to maximal such models predict the atmospheric mixing angle to be in 
order to determine whether the current or next generation of neutrino oscillation experiments is capable of favoring 
or disfavoring this particular approach to "explaining" the leptonic sector. Currently, the atmospheric angle is only 
loosely constrained by atmospheric neutrino data: —0.44 < cos20 2 3 < 0.44 at the three sigma confidence level 

(Eq. jnnu)).§ 

The current generation of long-baseline accelerator experiments is not expected to improve significantly on the 
current bound |23j , while studies of next-generation "off-axis" experiments have claimed that these are able to measure 
sin 2 2023 at the 1% level More detailed studies, however, are still absent [24|. Note that while a 1% measurement 
of sin 2 2023 sounds impressive, for a best-fit point sin 2 2# 23 = 1 it "only" translates into |cos2# 23 | < 0.1. Other 
constraint may be provided by future atmospheric neutrino studies |25l l2t| . or, ultimately, by neutrino factories, 
which may be able to measure cos 26> 2 3 at the few percent level [2^ • 

I have explored different neutrino Majorana mass-matrices that yield, at zeroth order, a v% state which is 0% v e> 
50% Vp, and 50% v r and are "perturbed" by a random Majorana mass-matrix. The magnitude of the perturbation 
is chosen such that one can also account for the solar mass-squared difference. In order to extract numerical results 



§ It is curious that the measurement of the atmospheric angle is significantly less precise than the measurement of the solar angle. 
Currently, 0.22 < cos26>i2 < 0.56 at the three sigma level. One of the reasons for this is the fact that strong solar matter effects render 
the survival probability of electron-type solar neutrinos proportional to sin 2 012 for a large range of solar neutrino energies, and hence 
the solar data is very sensitive to 0i2. 
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TABLE I: Order of magnitude for the average value of the leptonic mixing parameters for the different neutrino mass-textures. 
See Sec. 11111 for more detail. Also included are the expectations for an anarchical mixing matrix plllll. 



Texture (Sec.) 




| COS 26*23 | 


Solar Angle 


Normal Hierarchy lintV) 


^/Aro? 2 /Amf 3 


y/Am 2 2 /Am 2 3 


0(1) 


Inverted Hierarchy ijHIB) 


Am 2 12 /\Am 2 13 \ 


Am 2 12 /\Am 2 13 \ 


0(1) 


Inverted Hierarchy, Traceless ijHICl 


Am 12 /|Am? 3 | 


Ami 2 /\Am 2 13 \ 


sin 6>i2 ~ Arai2 / Am\ 3 


Inverted Hierarchy, Bi-maximallHIfD) 


Am? 2 /]Am? 3 | 


Am 2 12 /\Am 2 13 \ 


|cos26* 12 | ~ Am 2 12 /\Amj 3 \ 


Anarchy 


> 0.1 


0(1) 


0(1) 



from the random mixing matrices, I computed the absolute values of different observables, and replaced the order one 
numbers by their average values. This procedure allows one to easily define the magnitude of the perturbation and 
estimate how much | 1 (|t/ e 3|) is expected to deviate on average from 1/2 (0). A summary of the results obtained 
is presented in Table [I] 

For all mass-textures examined, the following features are observed. Both | ?7e3 1 and | cos 26*23 1 are, on average, of 
the same order of magnitude, as one should exp ect from a featureless perturbation (this is obviously not the case for 
all neutrino mass models. See, for example, [28] for a counter-example). Furthermore, they are of order Am 2 2 /| Ato 2 3 | 
for textures that yield an inverted mass-hierarchy or of order y Am 2 2 /| Am 2 3 | for textures that yield a normal mass- 
hierarchy. As far as the solar mixing angle is concerned, it is either "anarchical" (order one), in good agreement with 
the data (this is the case of mass-textures that yield v\ and Vi states that are mass-degenerate at zeroth order) , or 
doomed to be, on average, way too small/too large, in which case, the current approach safely does not yield a good 
fit to the data. 

What do these results mean? In the case of a normal mass-hierarchy, | cos 20231 is, on average, significantly away 
from maximal, such that proposed off-axis experiments should observe a small but significant deviation of sin 2 26*23 
from 1. Clearly, if they do not see such a deviation, it does not mean that mass textures that yield a normal mass 
hierarchy are "ruled out" - it simply means that they either "got lucky" or that they are not "perturbed by a random 
mass-matrix." 

In the case of an inverted mass-hierarchy, #23 is very close to maximal (and |t/ e 3| 2 is <SC 10 -2 , a rather unfortunate 
case for the off- axis experiments under consideration and reactor searches for |t/ e 3|) and no deviation from maximal 
should be observed in next-generation experiments. Note that the average upper bounds calculated above can be 
violated by a factor of a few due to random statistical fluctuations of the order one parameters, but are probably 
not violated by more than an order of magnitude even for more structured models. Keep in mind that the type of 
models under consideration here predict maximal mixing at zeroth order. Perturbations that are consistent with the 
solar data and that yield a much larger deviation of cos 26*23 from zero will violate this assumption. Therefore, the 
experimental determination that 100826*23! is not consistent with zero at the 10% level indicates that textures that 
yield maximal atmospheric mixing and an inverted mass-hierarchy are disfavored. 

Before concluding, I should mention that the analysis performed here is not, and is not meant to be, exhaustive. 
There are, potentially, other models that lead to maximal mixing in the atmospheric sector at zeroth order but are 
not well characterized by the mass-textures explored. I would also like to emphasize again that I only address textures 
that predict maximal atmospheric mixing at zeroth order. There are several models in the literature that "predict," 
say, large atmospheric mixing and zero U e 3 at zeroth order (one simple example is the so-called "semi-anarchical" 
texture of jl^j). As far as these are concerned, the expected deviation of the atmospheric angle from 7r/4 is order 1. It 
should be noted that, according to the results on Sec. lIII*A . it may not be easy to differentiate a texture like Eq. (|III.2|) 
from the semi-anarchical one if cos 26*23 is determined to be significantly away from - more detailed analyses would 
be required. 

In conclusion, neutrino oscillations provide our first and, currently, only evidence for physics beyond the Standard 
Model. The amount of information available regarding the nature of this new physics is, however, very small. As far 
as the neutrino sector is concerned, and in the case of neutrino oscillation experiments in particular, it is possible 
that the best one can accomplish is to measure precisely the elements of the leptonic mixing matrix, the value of the 
mass-squared differences, and the neutrino mass-hierarchy. Hence, we need to take maximal advantage of all these 
measurements in order to come up with the most complete and useful description of Nature. 

I tried to point out in this exercise that the deviation of the atmospheric mixing angle from maximal is, potentially, 
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an important observable, which may be able to test whether this maximal mixing phenomenon is somehow encoded 
in Nature via some fundamental principle, or whether it is simply accidental. It seems that even if maximal mixing 
is indeed fundamental at some energy scale, it is still true that, in the case of a normal mass hierarchy, it is possible 
that next-generation oscillation experiments will see a deviation of cos 2#23 from 0, but more quantitative analysis 
of these experimental proposals are required in order to better quantify this statement. On the other hand, in the 
case of and inverted mass-hierarchy, cos 2#23 = seems to be a robust prediction for any next-generation neutrino 
oscillation experiment, given that these probably cannot measure cos 2#23 at better than 10%. 

NOTE ADDED: After the completion of this manuscript, |30j became publicly available at the preprint archive. 
It contains a discussion of the effects of a "perturbation" on a zeroth-order bi-maximal leptonic mixing matrix (the 
presentation and motivation are different from the one here. See |30f for details), and discusses relations between the 
atmospheric mixing angle and \U e3 \. It also contains some of the results presented in Sec. HHP . 
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APPENDIX A: COMBINING ORDER ONE RANDOM ENTRIES 



In order to obtain numerical estimates for the products and sums of the random order one parameters a, 6, c, d, e, / 
defined in Eqs. i|II.2IIII.l|) . I follow the prescription below: 

i- Assume that a, b, c, d, e, / are uncorrelated random variables, with a Gaussian probability density function centered 
at zero and with identical widths, 



V27 



■. exp 



2a 2 



x = a, b, c, d, e, /. 



(A.l) 



The probability that x is between xi and x u (xi > x u ) is denned to be f*" g{x)dx. The width a is defined such that 
the average absolute value of the trace of Mi = mi. In the 2x2 case, Eq. (|II.2() . a — y/n/2 (see below), while in the 
3x3 case, Eq. iftlTTjl . a = 

ii-Replace combinations of order one parameters by their average values according to their probability density 
functions. For example 



\a + b\ = / \a + b\g(a)---g(f)da---df = / \a + b\g(a)g(b)dadb, 



2ira 2 



27TCT 2 



27TCT 2 



b | exp 



exp 



-±((a + b) 2 + (a-b) 2 ) 



d{a + b)d{a-b) 



x 

Aa 2 



2 \ ^ 

2 ' 



-4a' = 



2a 



(A.2) 
(A.3) 
(A.4) 
(A.5) 



In the 2x2, \a+ b\ = 1 (by definition) while in the 3x3 case, \a + b\ = \/2/3. 

For other simple averages of sums or products, one can follow a similar procedure and obtain analytic results quickly. 
One other example is 



b - 2c) 



1 



(27T ( J 2 ) 3 /2 

2V67RT 2 
(27TCT 2 ) 3 / 2 



(a + b - 2c) 2 exp 



(\( a + b - 2c ) 2 + \( a - b f + l( a + b + c ) 2 



x 2 exp 



^2ct 2 V6 
2 ^ dx 



\f%ita 2 



12a 2 J 6 3(2na 2 ) 3 / 2 



(6a 2 ) 3 / 2 V2^, 



= 6a'. 



dadMc(A.6) 

(A.7) 
(A.8) 



Hence, in the 2x2 case, (a + b — 2c) 2 = 3tt/2, as quoted in Sec. ITT1 

Another, related, possibility would be not to compute average values, but instead identify "confidence level ranges" 
for the various outputs as a function of the original width of the individual probability distributions, a. For example, 
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Eq. (|II.6(1 could be reinterpreted as: the probability distribution for Am 2 /mQ is centered at zero and has width \/2<r. 
After defining the value of a (by requiring, say, that Tr(Mi) < mi at the one-sigma confidence level), by requiring that 
the experimentally observed value is within some pre-defined confidence level, one can extract the value of mi/mo, as 
was done for the average-value procedure outlined above. Note that the requirement that \a + b\ = 1 translates into 
a width o a +b = Hence, the probability that a + b is between —3 and 3 is around 77%. 
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